Finite hadronization time and unitarity in quark recombination model 
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The effect of finite hadronization time is considered in the recombination model, and it is shown 
that the hadron multiplicity turns out to be proportional to the initial quark number and unitarity 
is conserved in the model. The baryon to meson ratio increases rapidly with the initial quark density 
due to competition among different channels. 
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In last few years, the quark recombination (or coales- 
cence) model 0, 0, IE has aroused intense interest in 
the communities of both theorists and experimentalists 
in the field of ultra-relativistic heavy ion collision physics, 
because of its ability to explain novel phenomena, such 
as unexpected high proton over pion ratio at px about 
3 GeV/c @ and the constituent quark number scaling 
of the elliptic flow 6], that are observed in RHIC ex- 
periments. Although different physics considerations are 
input in the investigation of final state particle spectra 
in different implementations of the quark recombination 
model by different groups, there are a lot of common 
points |3- All implementation of the model, except the 
latest attempt 0, considered only quark and antiquark 
degrees of freedom and gluons are assumed having been 
converted into quarks and antiquarks before hadroniza- 
tion. Every group assumed implicitly that hadronization 
occurs instantly. In heavy ion collision processes there 
are plenty of soft thermalized partons and some hard par- 
tons which will evolve into showers of semihard partons 
0. Because of the absence of gluon degree of freedom 
in the model, quarks (and antiquarks) at the moment 
just before hadronization must have been dressed up as 
constituent quarks ready for recombining into final state 
hadrons. It has been shown that the recombination of 
thermal quarks dominates in Au+Au collisions at RHIC 
energies for low pt pion production and the traditional 
fragmentation, which can be interpreted as a recombina- 
tion process from shower partons originated from a hard 
parton 0, dominates at very high pt, while the thermal- 
shower recombination is important for moderately high 
pt. 3GeV/c < pr < 8GeV/c in central Au+Au collisions 
101. In this way, all hadron production can be consid- 
ered consistently in the framework of quark recombina- 
tion model. In every implementation of the quark recom- 
bination model, the transverse momentum spectrum of 
mesons at mid-rapidity can be written, after some alge- 
bra, as 
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where a factor S(pt — Pi — p-i) is included in the recom- 
bination function R M (pi,P2,Pt) to ensure momentum 
conservation and F{p\,p2) is associated with the joint 



quark-antiquark momentum distribution. When we are 
interested in the total multiplicity of a kind of meson, we 
can consider the contribution from thermal-thermal re- 
combination only, because most produced hadrons are in 
the low transverse momentum region where pure thermal 
recombination dominates. In this range of transverse mo- 
mentum, the joint distribution for quark-antiquark pair 
can be written as F(pi,p 2 ) — V p 2 ' fi{pi)f2{p2) with V 
the spatial volume of the partonic system and p the ther- 
mal parton density just before hadronization when the 
(anti)quark transverse momentum distributions /i,2 are 
normalized to 1. A simple conclusion from Eq. Q is 
that the total yield of the meson is proportional to the 
square of quark density just before hadronization. Simi- 
larly, one can conclude that the yield of a kind of baryon 
should be proportional to the cubic of the quark den- 
sity. On the one hand, if one considers production of 
pentaquarks along the same line, one has to say that its 
yield is proportional to p 5 and higher in heavy ion colli- 
sions than in elementary ones. This conclusion seems to 
contradict to current experimental observations. On the 
other hand, that the yield in a given coalescence channel 
scales quadratically/cubically with the constituent num- 
ber violates unitarity, as argued in 00,0, because the 
total number of hadrons should be linearly proportional 
to that of constituent quarks just before hadronization. 
Apparently, the total number of constituent quarks in fi- 
nal state hadrons should be equal to the quark number 
just before hadronization. Considering above, one has 
to conclude that something must be changed in current 
formulism of the recombination model. 

An important observation in this paper is that the 
numbers of constituent quarks are decreasing at finite 
rate during the hadronization process if one assumes that 
the hadronization lasts a nonzero period of time. Then 
one should study the effect of finite hadronization time 
on the yields of hadrons in the recombination model. For 
this purpose, one can interpret the recombination func- 
tion R M as the average production rate of meson from a 
given pair of quark and antiquark and write the joint two- 
quark distribution as a function of time explicitly. Then 
the right hand side of Eq. £Q) should be dN M /p T dpTdt. 
Here time t can be defined in the laboratory frame. In 
this paper, we will show that the introduction of the fi- 
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nite hadronization time in the recombination model can 
ensure unitarity. 

Instead of considering the spectra of different hadrons 
in the framework of quark recombination models, we only 
investigate the yields of mesons and baryons, as done in 
earlier work of the coalescence model |l3|. For simplic- 
ity, we assume in this paper that quarks and antiquarks 
have the same density in the almost baryon free central 
rapidity region and that the shapes of the all transverse 
momentum distributions involved remain unchanged in 
the hadronization process. The strange contents are not 
considered in this paper. If the phase transition from par- 
tons to hadrons is of first-order, there is a boundary and 
hadrons are produced only on the surface while the par- 
tons inside have a constant temperature and/or density. 
If, on the other hand, the phase transition is of second- 
order, no boundary between parton state and hadron 
state can be defined, and hadronization takes place in 
the whole system at the same time. Most event genera- 
tors for high energy collisions adopt the second scenario 
in the codes. Such a scenario was also taken in [J^j. In 
ultra-relativistic heavy ion collisions the produced par- 
tonic system expands. Then the phase transition may 
be neither first-order nor second-order. If the expansion 
retains to the last stage of the evolution and the phase 
transition is likely first-order, there are two competing 
trends on the change of the volume. The expansion will 
make the partonic volume larger, and the hadronization 
tries to shrink the system. Thus the hadronization dy- 
namics should, generally speaking, be connected with the 
hydrodynamical calculations. This cannot be included in 
this paper. Now we simply consider two limiting cases 
for the evolution of the volume in hadronization process. 

The first case we consider is for fixed volume for par- 
tonic system. In other words, the expansion of the sys- 
tem is assumed to be compensated by the shrinking from 
hadronization. Because the hadronization time is very 
short, this case can be regarded as the limit of extremely 
rapid expansion. 

First of all, we assume that only one species of meson 
can be produced in the process. Then one can integrate 
over the momentum of produced hadrons and get the 
production rate of the meson as 
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and the information of both the shape of quark distribu- 
tion and the recombination function is encoded in A. No 
reverse term for the process from hadron to quarks is con- 
sidered here, because (i) most of the produced hadrons 
will move quickly to the detectors, thus the hadron den- 
sity in the reaction zone is much smaller than that for 
quarks, and (ii) at the energy scale corresponding to the 
phase transition (about 170 MeV), it is very unlikely for 
hadrons to dissolve into quarks and antiquarks in the 
interactions. As a consequence, the interactions among 



the produced hadrons will be elastic and no lose term 
appear in Eq. (2). With the production of the meson, 
the constituent quark number decreases accordingly, and 
we have 
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dt 
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This equation comes from the model assumption that 
just before hadronization all partons have been converted 
into constituent quarks ready for recombination into fi- 
nal state hadrons and no more constituent quarks can be 
created. Finally, all those N q + Nq constituent quarks 
have to be inside some hadrons, because of color confine- 
ment. Thus as time goes to infinity, p must go to zero 
from initial density po = N q /V just before hadroniza- 
tion. Equations corresponding to Eqs. (2) and (3) can 
be found for any chemical reaction A + B — > C which 
is obviously a combination process similar to the process 
q + q — > 7r studied in the quark recombination model. 
For chemical reactions the amount of final product al- 
ways depends linearly on the initial density. From above 
two equations, one can easily deduce the yield as 
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V Po = N q 



(4) 



Therefore, the total multiplicity of the meson turns out to 
be equal to the initial number of constituent quark, and 
unitarity is maintained in the recombination process. A 
careful reader may have noticed that the parameter A 
introduced in Eq. (2) is not relevant to the total mul- 
tiplicity of the meson, as shown in Eq. (JIJ. This is 
expected from quark number counting. 

If only one species of hadrons consisting of n con- 
stituent quarks is considered, the total multiplicity N n 
can be shown, along the same line as for mesons, to be 



N n 



2V Po 2N q 



(5) 



Again the total multiplicity is proportional to the initial 
number of constituent quark and there is no violation of 
unitarity. 

In real high energy collisions, many different species of 
particles are produced. In the recombination model, the 
production rates for different hadrons depend on quark 
distribution differently. For mesons the production rate 
depends on quark density quadratically, while for baryons 
the dependence is of cubic. So there exist competitions 
among different channels and the effect of competition 
depends on the initial number of quarks. To get a pa- 
rameter free expression, we assume in this investigation 
that only non-strange mesons and (anti)baryons are pro- 
duced. Then we have 
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= A 2 V P 3 (t) , 
fdN M 
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(6) 
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with initial conditions N M > B = at t = 0, p(0) = p . 
There are exactly the same number N B anti-baryons. 
The third expression in last equation ensures the unitar- 
ity in hadronization process. Because unitarity condi- 
tion was not taken into account in former applications 
of the recombination model, the result in this paper will 
not be reduced to the old ones even in the limit of zero 
hadronization time. From the reaction-rate theory 
A\ and A 2 are proportional to the corresponding cross- 
sections and the degeneracy factors. Above equations can 
be solved numerically. First we can define dimensionless 
quantities u — 3A 2 po/ A\,r = p(t)/po,T = Aipot. Here 
variable u has encoded initial quark density and the com- 
petition between meson and baryon channels. With these 
scaled variables above equations can be normalized to 
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(7) 
(8) 



where the scaled quark density r(u, r) at scaled time r is 
the solution of the equation 



(9) 
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FIG. 1: Behavior of scaled density as a function of scaled time 
for fixed u = 1 and 5. 

The relative density r are shown in Fig. 1 as a function 
of scaled time r for u = 1 and 5, respectively. For larger 
u, r decreases more quickly with t. This is not surprising 
because of the quadratical and cubical dependence on r 
of the rates for meson and baryon production. From the 
fact that the constituent quark density decreases with 
r very rapidly, we can see that the average production 
time for baryons is shorter than that for mesons. Or in 
other words, baryons are produced earlier than mesons 
in the hadronization process. It has been shown that the 



pion size measured by interferometry increases with the 
pion multiplicity [la, steadily with bombarding energies 
in similar systems |l7j |. For protons, the measurements 
of NA49 [13 together with measurement at lower ener- 
gies suggested a very weak dependence of the two-proton 
correlation function on bombarding energies. If protons 
and pions share the same collective expansion, smaller 
source size means earlier production. Our results agree 
qualitatively with those experimental conclusions. 




FIG. 2: Dependence of the scaled multiplicities of mesons and 
baryons on the scaled initial density u. 

Once we know the relative quark density as a function 
of scaled time r, one can integrate Eqs. (0 and JHJ to 
get the multiplicities of mesons and baryons produced in 
the hadronization process. There exists an unknown con- 
stant factor VAi/A 2 . In the following, we define scaled 
multiplicities of mesons and baryons as 



N M,B 



ZN M ' B A 2 



(10) 



uAiV 

One can easily see that N s is proportional to hadron yield 
divided by the quark density just before hadronization. 
N^ 1,3 are shown in Fig. 2 as functions of u. As the initial 
quark density increases, the scaled yield for mesons de- 
creases, indicating that the meson yield increases slower 
than the constituent quark density does. Meanwhile the 
scaled yield for baryons increase quickly at small u, there- 
fore the baryon yield increases faster than the quark den- 
sity. The difference between the trends of Nm and Nb as 
u increases can be shown more clearly by the baryon to 
meson ratio N B /N M , as in Fig. 3. The ratio increases 
with u almost linearly in the range shown. In the near 
future at LHC, it can be predicted from Fig. 3 that the 
baryon to meson ratio will become even larger, since the 
initial constituent quark density (or u) will be higher. 
If the scaled constituent quark density can become high 
enough (about u = 8), the yield of baryons will be about 
the same as that for mesons, and the baryon to meson 
ratio can be more than 2 at u — 20, if it can be achieved. 
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FIG. 3: Baryon to meson ratio N B /N M as a function of scaled 
density u. 



For the second case we assume that there is no expan- 
sion at all. Thus the system shrinks to from initial 
volume Vq while the parton density remains the same in 
hadronization. Since p is now a constant, we get from 



Eq. © 



(11) 



In this case the hadronization is much faster than in the 
first case, because the volume of parton phase decreases 
exponentially. The multiplicities of produced mesons and 
baryons are 



N M = 



A\N„ 



A ± + 3A 2 p 



N B = 



A 2 N qP 
Ax + 3A 2 p 



(12) 



Then we have the baryon to meson ratio N B /N M = u/3. 
Now the ratio turns out to be proportional to the initial 
parton density and is larger than the corresponding val- 
ues for the first shown in Fig. 3. 

In real ultra-relativistic heavy ion collisions the parton 
system is always expanding in hadronization. The expan- 
sion will slow down the shrinking of the parton volume 
in hadronization. Thus the two cases considered in this 
paper are the two limiting evolutions of the system. The 
real dependence of hadron yield will be in between the 
results for the two limiting cases. 

One may have noticed that we only considered in this 
paper the contribution to particle production from the 
recombination of plenty soft quarks, which is the most 
important at low pr region. This discussion may be ap- 
plied to particle production at SPS, RHIC and LHC en- 
ergies with different initial parton density or volume as 
long as a dense parton phase is created. Thus the rapid 
increase of baryon to meson ratio can be a signal of the 
formation of the dense medium of soft partons, and the 
density of the medium can be deduced approximately 
from the measured p/ir ratio. The initial parton density 



just before hadronization may depend only on the collid- 
ing energy. So the baryon to meson ratio for collisions at 
fixed center of mass energy will be proportional to the ra- 
tio of the reaction rates, or in other words, proportional 
to the ratio of the degeneracy factors, as claimed by the 
statistical hadronization model [l9|. 

It should be pointed out that the shape of constituent 
quark distributions can, in reality, be changing in the 
hadronization process. As a result, above results may 
be modified. However, one can easily imagine that the 
shape changes more slowly than the normalization of the 
distribution. Then the modification will not be too big. 
In Eqs. (7) and (8) the shape dependence can be canceled 
partially in Ax/A 2 . So, qualitatively, the results from this 
investigation will not be changed when the evolution of 
shapes of quark distributions is taken into account. 

In conclusion, we showed in this paper the impor- 
tance of the evolution of constituent quark density in the 
hadronization process. When the evolution is taken into 
account, unitarity can be maintained, and the baryon to 
meson ratio is shown to increase monotonically with the 
initial constituent quark density. 
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